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Abstract 



The main topic of this paper is various "hyperbolic" generalizations of the Edmonds- 
Rado theorem on the rank of intersection of two matroids. We prove several results in this 
direction and pose a few questions. We also give generalizations of the Obreschkoff theorem 
and recent results of J. Borcea and B. Shapiro. 

1 Introduction, Spectral inequalities 

Consider a homogeneous polynomial p(x\, x m ) of degree n in m real variables. Such a 
p is called hyperbolic in the direction e (or e- hyperbolic) if for any x E R m the polyno- 
mial p(x — Ae) in the one variable A has exactly n real roots counting their multiplicities. 
We will assume below that p(e) = 1. Denote an ordered vector of roots of p(x — Ae) as 
X(x) = (\i(x) > \2(x) > ...A n (x)). It is well known that the product of roots is equal to p(x). 
Call x E R m e-positive (e-nonnegative) if X n (x) > (A n (x) > 0). The fundamental result ^H] 
in the theory of hyperbolic polynomials states that the set of e-nonnegative vectors is a closed 
convex cone. A fc-tuple of vectors (x\, ...Xk) is called e-positive (e-nonnegative) if Xi, 1 < % < k 
are e-positive (e-nonnegative). 

Below, we denote the closed convex cone of e-nonnegative vectors as A^, and the open 
convex cone of e-positive vectors as C e . It has been shown in ^B] (see also j2HJ) that an e- 
hyperbolic polynomial p is also d- hyperbolic for all e-positive vectors d E C e . 

Let us fix n real vectors X{ E R m , 1 < i < n and define the following homogeneous polyno- 
mial: 



Equivalently, the p-mixed value M p {x\, ..,x n ) can be defined by the polarization (see [22]) : 




(1) 



Ki<n 




(2) 






(3) 



6ie{-l,+l},l<i<n l<i<n 



Ki<n 



Let us denote as Ik tn the set of vectors r = 
and J2i<i<k r i = n. 



(ri,...,rk) with nonnegative integer components 



Let us associate with any such vector r an n-tuple of m-dimensional vectors X r of n copies of 
Xi(l < i < k). Notice that 

X r = (yi, ...,y n );y% G {x\, ...,x k }, 1 < i < k. 

It follows, for instance from the polarization identity (3), that 

P Xl> .., Xn (a 1 ,...,a n )= ]T II «?M p (X r ) 1 (4) 

r£/„,„ Ki<n 1 ll<i<n r »- 

For nonnegative tuple X = (xi, ..,x n ), define its capacity as: 

Cap(X) = inf ^...^(ai, at n ) (5) 

a « >0 'rii< l <„ a «= 1 



Example 1.1: Probably the best known example of a hyperbolic polynomial is 

P{a ,...,a k ) = Det{ £ o^) (6) 

0<?<fc 

where Ai , < i < k are hermitian matrices and the linear space spanned by Ai , < i < /c 
contains a strictly positive definite matrix: J2o<i<k Pi^i = B >~ 0. This polynomial is hyperbolic 
in the direction (3 = ((3i, (3 k ). We can assume wlog that B = I and that /3 = (1, 0, 0, 0). 
In other words, after a nonsingular linear change of variables 

P(a , ...,a k ) = Det( a i B i) ( 7 ) 

0<i<fc 

where the matrices Bi,l < i < k are hermitian and Bq = I. 

Let A\...A n be n x n matrices. 
• x n A n ) (8) 



Whether or not the cone N e of e-nonnegative vectors allows a semidefinite representation is a 
major open question in the area. In the case of three variables this question was recently settled 
in [21] which is a rather direct application of [Q. In this three variables case if e = (1, 0, 0) then 
any e- hyperbolic polynomial has a determinantal representation (7), in fact, it even has one 
with real symmetric matrices B{. 

Proposition 1.2: Consider a homogeneous polynomial p(x%, x m ) of degree n in m real 
variables which is hyperbolic in the direction e. For any pair of m-dimensional vectors x, y G R m 
there exist a pair of n x n real symmetric matrices A, B such that \{ax + by) = X(aA + bB), 
where a,b £ R and A(M) is the ordered vector of eigenvalues of the matrix M. 



In this case mixed forms are called mixed discriminants. 
The number n 

D(Ax...A n ) = — det(siAi + • 

OX\...OX n 

is called the mixed discriminant of A\...A n . | 
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Proof: Consider the following hyperbolic in the direction (1,0,0) polynomial Q{x\, X2, X3) = 
p{x\e + X2X + x^y). Then there exists two n x n real symmetric matrices A, B such that 
p{x\e + X2X + x^y) = Det(x\I + X2A + X3B). It follows that for real a, b the roots of ax + by 
coincide with eigenvalues of aA + bB and our proof follows directly. | 

Propositin 1.2 allows to "transfer" many "spectral" statements, known for real symmetric 
matrices, to the context of general hyperbolic polynomials. (We assume that conditions of 
Proposition 1.2 are also satisfied in Corollary 1.3 below). 

Corollary 1.3: 

1. Consider a symmetric (i.e. invariant respect to all permutations of variables) f(yi, ..,y n ) : 
X — » R , where X C R n and either X = R n , either X = R™ (nonnegative orthant) or 
X = R+ + (positive orthant). Define f(x) = f(Xi(x),X2(x)...X n (x)) where either x £ R n , 
either x £ N e or x £ C e correspondingly. If f is convex on X then f(cc) is also convex on 
either R n , either N e or C e correspondingly. (Most recent "hyperbolic inequalities" papers 
(131) Cffi etc.) are simple corollaries of this statement.) 

2. Recall that for two nx n hermitian A,B there is a complete polytope description (Horn's 
inequalities) of all possible triplets of vectors (\(A + B) , X(A) , X(B)) such that: 

52HA + B) <5>i(A) + J2 \k(B), 

where T,S,U are some suitable subsets of {1, ...,n} of the same size. 

We get from Proposition 1.2 that for any two vectors x, y £ R m the ordered vectors of roots 
(X(x + y), X(x), X(y)) satisfy all Horn's inequalities. In particular they satisfy the Lidskii 
property : there exists a doubly stochastic nxn matrix D such that X(x+y) — X(x) = DX(y). 
(This settles Open Problem 3.6 posed in \l0l). 

Remark 1.4: In the very same way one can obtain results on self-concordance. Indeed, results 
on self-concordance are of the following nature: consider some symmetric function /, and for 
a pair x,y £ R m define F(t) = f(x + th). Assume that x belongs to some cone, usually it is 
a cone of positive vectors |llj . One needs that F is convex and satisfies some inequalities for 
derivatives of F at zero: 

| F (3)| < 2(F (2) ) 1.5. | F (D| <^ (jF (2)).5 

Again, if these inequalities and convexity hold for hermitian nxn matrices, then we get the 
same stuff for vectors/e-positive vectors/ e-positive vectors satisfying p(x) > a > using the 
hyperbolic polynomial p{x\e + X2X + x^h) = Det{x\I + X2A + x%B). 

The class of inequalities which follow from Proposition 1.2 is larger then the class of sym- 
metric convex inequalities from |10j . For a complex matrix A, consider its singular values 
01 > a 2 > ... > a n > 0. Define f{A) = i/E^- Then f(A + B) < f(A) + f(B) HE], |2H- 
Other inequalities of this type, which are obtained using optimal nonholonomic control, can be 
found in j^H]. And all of them are non-convex, all of them can be "transfered" to hyperbolic 
polynomials. 
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Many other things become much more transparent using the polynomial (in three real variables) 
p(x\e + x 2 x + X3/1). For instance, the mentioned above fact that e- hyperbolic polynomial p is 
also d- hyperbolic for all e-positive vectors d £ C e and the equalities 

N e = N d ,C e = C d 

We will get more applications of this "trick" (i.e. using hyperbolic polynomial p(xie+x 2 x+x 3 y) 
in three variables) in Section 3. I 

2 Inequalities for mixed forms, Combinatorics of mixed forms 

In this section we will try to understand if another important class of inequalities, which is valid 
for determinantal polynomials (7) , can be "transfered" to general hyperbolic polynomials. We 
recall below some facts about mixed discriminants, which are mixed forms corresponding to the 
determinantal polynomials. 

Alexandrov-Fenchel inequalities. Consider an n-tuple of positive-semidefinite n x n her- 
mitian matrices (A\, ..,A n ). Then the mixed discriminant satisfies the following 
(hyperbolic) inequality : 

D(A 1 ,A 2 ,A 3 , ..,A n ) > ^D(A 1 ,A 1 ,A 3 ,..,A n )D(A 1 ,A 2 ,A 3 ,..,A n ) (9) 

This inequality holds also for mixed forms M v {x\, ..,x n ), where p is e-hyperbolic polyno- 
mial of degree n , and (xi, ..,£„) are e-nonnegative vectors [25] . 

van der Waerden inequality Call an n-tuple of positive-semidefinite n x n hermitian ma- 
trices (Ai, ..,A n ) doubly stochastic if 

tr(Ai) = 1, 1 < i < n; ^ Aj = I (10) 

l<j<n 

Then 

111 1 r?' 

D(A 1 ,A 2 , A 3 , .., An) > D(-I, -I, -I, —I) = — (11) 

n n n n n n 

Moreover the inequality is strict if (Ai,A 2 ,A 3 , ..,A n ) / (^J, ^1, ±1, ±1). 

The inequality (11) was conjectured in |7j and was proved for the real case in [20] : the 

hermitian case and uniqueness were proved in |17j . 

Scaling Here is the version of (11) which does not require doubly stochasticity. Similarly to 
(5) define 

Cap(Ai,A 2 ,A 3 ,..,A n ) = inf Det{a x Ax, a n A n ) (12) 

^ >o -ni< l <„^= i 

Then 

n n 

D{A 1 ,A 2 ,A 3 ,..,A n ) < Cap(A 1 ,A 2 ,A 3 ,..,A n ) < — D(A 1: A 2 , A 3 , .., A n ) (13) 

n! 
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The proof of this inequality in |2L)| . |21| requires a matrix scaling. Left inequality in (13) 
holds also for mixed forms of nonnegative vectors for general hyperbolic polynomials. 

"Concavity" of ln(Cap) We present below a general result, i.e. which holds for general 
hyperbolic polynomials. 

Lemma 2.1: Consider an e-nonnegative tuple X = (x±, ..,x n ). For a vector r = 
(ri,...,r n ) £ I(n,n) define f(r) = ln(Cap(X r ). The function f(r) is concave on I(n,n). 
I.e. ifr Q = J2i<i<k a i r i> where J2i<i<k a i = l;Oj>0,l<i<A; and n £ l(n,n),0 < i < 
k, then the following inequality holds : 

Cap(X ro )> J] Cap(X r J a * (14) 

l<i<k 

Proof: We will use a known recent result ^U] that the functional ln(p(x)) is 

concave on a positive cone C e (see also Proposition 1.2 and Section 4 ). Fix a vector 
(ai, a n ) : a>i> 0, Ui<i< n a i = 1 - 

First, let us consider some z = (zi,...,z n ) £ I(n,n). An easy application of the arith- 
metic/geometric mean inequality gives that 

Cap(X z ) >dittp(Yl > d IT a ? ( 15 ) 

l<j<n l<i<n 

Let T{ = (rj i,ri 2j •••> r i,n)jO < i <k. It follows from (15) that 

ln(p( y~] ctjrijXj)) > ln(ctj)rij + ln(Cap(X ri )), 1 < « < 

Multiplying the ith inequality by etj and adding afterward we get that 

ailn{p( Y a j r i,j x j)) > ln(aj)r j + ^ ailn(Cap(X n ) 

l<i<k l<i<" l<7<n l<i<fc 

Using the concavity of ln(p(.)) and (15) we finally get that 

Cap(X ro )> J] Cap(X ri ) a ' 
i<j<fc 

I 

Edmonds-Rado theorem and Newton polytopes The following result is a direct corol- 
lary of the famous Edmonds-Rado theorem on the rank of intersection of two matroids 

M ■■ 

Consider an n-tuple of positive-semidefinite n x n hermitian matrices (Ai, ..,A n ). Then 
the the mixed discriminant D(Ax, A%, A3, .., A n ) > iff Rank(^ i&s Ai) > \S\ for all 
S C{l,2,...,n}. 
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Consider a polynomial Q( a i, a n ) = (Det{oi\Ai + ... + a n A n ). Let us denote as Ik, n the 
set of vectors r = (n, with nonnegative integer components and X)i<i<fc r i = n - 

Associate with an integer vector r € /(n, n) an n-tuple of matrices A r consisting of j*j 
copies of Aj(l < i < n) and denote by D(A r ) the corresponding mixed discriminant. 
Then 

Q(a 1 ,...,o n )= £ 1] ^£>(X r ) 1 (16) 

r€l n ,„ l<i<n lll<i<n r «- 

Therefore the support of Q, supp(Q) = {r £ I n ^ n : D(K r ) > 0}. 
It follows from Edmonds- Rado theorem that 

CO(supp(Q)) n I(n, n) = supp(Q), (17) 

where CO{supp(Q)) is a convex hull of supp(Q) , i.e. the Newton polytope of the poly- 
nomial Q . 

The inequality (13) and Lemma (2.1) give a more precise statement : If ro = J2i<i<k a i r ii 
where J2i<i<k a i = 1) > 0, 1 < i < k and G I(n,n), < i < k, then the following 
inequality holds: 

D(X ro ) > J] D(X n r^ (18) 

It is interesting to notice that proofs of as (11) as well (13) do not use Edmonds-Rado 
theorem. 

Now we are ready to ask relevant questions for hyperbolic polynomials. 

Question 1. Consider a homogeneous polynomial p(yi, ...,y m ) of degree n in m real variables 
which is hyperbolic in the direction e and the corresponding polynomial in n real variables 

P Xl ,..,x n (ai,-.,a n ) =p( a i x i) 

l<i<n 

where xi, -.,x n G R m are e-nonnegative. Is it true that 

CO(supp(P xl! .. ;Xn )) n/(n,n) = supp(Q)7 
(Recall that the convex hull CO(supp{P xlt .^ Xn )) is the Newton polytope of the polynomial 

P X l,..,x n -) 

If the answer is "yes" then we get an analog of Edmonds-Rado theorem for hyperbolic 
polynomials. To state this, define the p-rank of x £ R m as Rank{x) = \{i : Xi(x) ^ 0}|. 
Then, the statement is that for e-nonnegative tuples X = (xi, ..,x n ) the p-mixed value 
M p (xi, ..,x n ) > iff Cap(X) > 0; or equivalently, iff Rank{^ i&s Xi) > \S\ for all S C 
{l,2,...,n}. 

This (conditional) result follows from the following Proposition 2.2 and Proposition 
2.3 . 

Proposition 2.2: Consider a homogeneous e-hyperbolic polynomial p(.) of degree n in 
m real variables. Let X = (xi,..,x n ) be an e-nonnegative tuple. Then Cap(X) > iff the 
following generalized Edmonds-Rado condition holds: 
Rank(J2ies x i) — \&\ f or a ^ & c {!> ^ ; —,n}. 
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Proof: We will use two known facts: 

Fact 1. \k(x + y) > \k(%) provided y is e-nonnegative; if x is e-nonnegative and the scalar 
a > then Afc(ax) = aXk(x), 1 < k < n. 

Fact 2. p{x) = p(e) Ili<fc<n ^k(x); if x, y are e-nonnegative then p{x + y) > p(x). 



If the generalized Edmonds-Rado condition holds then for any subset S of cardinality k 
we have the inequality 

i&S 

Using Fact 1 and the same argument as in the proof of Theorem 1.14 from Jl], we get 
that for e-positive (zi, 1 < i < n) the following inequality holds: 

p{ z % x i) > II rkZk 

l<i<n l<k<n 

In other words, Cap(K) > Y\\<k< n r k- 

Suppose that the generalized Edmonds-Rado condition does not hold ; or, wlog , (5Zi<i<fc x i 
0. Also, assume wlog that all vectors e — X{ are e-nonnegative , where e is a vector (direc- 
tion) used in the definition of hyperbolicity (recall that p(e) = 1). Choose = a > for 
1 < i < k and zi = b > for k + 1 < i < n. Using Fact 2, we get the following inequality: 

p{ Y z i x i) - P( X! aXi + ( n ~ k)be) 

l<i<n l<i<k 

By our assumption, an e-nonnegative vector J2i<i< n x i nas a ^ m °st k — 1 nonzero roots: 
f > r± > ... > r^_i > 0. Therefore 

p{ J2 ax i + ( n ~ k ) be ) = b n ' k+1 II an + b< b n - k+1 (ar + b) k ~ l 

Ki<k Ki<k-1 



Finally, 



P(J2i<i<n z i x i) {ar + b) k x b 



]ll<j<n z i ak 

For a fixed a > the right side of the last inequality converges to zero if (positive) b 
converges to zero. Therefore, Cap(X.) = 0. I 

Proposition 2.3: For any vector r = (ri,...,r n ) £ I(n,n) the capacity Cap(X r ) > iff 
r G CO{supp{P xl _ Xn )). 

Proof: It is an easy application of convexity of the logarithm on the positive orthant 
and the Hahn-Banach separation theorem; all what is important is that the coefficients 
of P X i,..,x„ are nonnegative. I 

Definition 2.4: Consider a polynomial p(y\, ...,y m ) of degree n in m real variables hy- 
perbolic in direction e and assume that P{e) = 1. Call an n -tuple X = (xx,..,x n ) 
of m-dimensional real vectors (i-doubly stochastic if Xj, 1 < i < n are e-nonnegative, 
J2i<i<n x i = d £ C e and tr^Xi) = 1, 1 < i < n , where tr^x) is the sum of roots of x in 
the direction d £ C e , i.e. roots of the equation p(x — td) = 0. I 
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Proposition 2.5: If an n-tuple X is d-doubly stochastic then Cap(X.) = p(d) 

Proof: Since p(J2i<i< n x i) = P(d) hence Cap(X.) < p(d). It remains to prove that 
p(J2i<i<n eaix i) — p(d) if J2i<i<n a i = an d 0-i are real. The functional g(a\, a n ) = 
p(J2i<i<n eaix i) i s convex on R n (even log(g(a\, a n )) is convex) as all coefficients of 
Pxi,..,x„ are nonnegative. 

Thus we need to check, similarly to the proof of Lemma 3.6 in j2J> that the gradient of 
g evaluated at the zero vector (0, 0) is a constant multiple of (1, 1). But at the zero 
vector, 

d d 
-—(g(a 1 ,...,a n )) = —p( V (1 + aijxi), 1 < k < n 
oab oai . ~ 

K 1 l<t<n 

Using d-double stochasticity, we get that 

p{ 0- + a i) x i) = P(d + aiXi ) 

l<i<n l<i<n 

Finally, at the zero vector, we get that 

d d 
■K—{g(ai,...,a n )) = -^—p{d + a k x k ) = p(d)~ 1 tr d {x k ) =p(d)~ ,1 <k<n 
oa k da k 

Therefore, the zero vector (0, 0) is a global (not always unique) minimum of g(a%, a n ) 
on the hyperplane (oi, ...,a n ) : J2i<i<n a i = 0- This means that 

p( Y OLiXi)>p(d) Y[ ot^oa > 0, 1 < i < n 

l<i<n l<i<n 

Thus Cap(X) = p(d). | 



Remark 2.6: Perhaps, Proposition 2.5 can be strengthened to the following statement: 
let A = (Ai,...,A n ) be roots of J2i<i< n c i x i m t ne direction d £ C e , where the tuple 
(xi, ...,x n ) is d-doubly stochastic and (q, 1 < i < n) are real numbers. Then there exists 
a doubly stochastic n x n matrix D such that A = CD, where C = (c\, c n ). 
I 



Question 2 Define the van der Waerden constant of a hyperbolic polynomial p(yi, ...,y m ) of 
degree n in m real variables as 



VDW(p) = inf 



M p (xi, ..,x n ) 
Cap(xi, ..,£„) 



where the infimum is taken over the set of tuples (xi,..,x n ) of e-positive vectors. It is 
easy to see that VDW{p) < Is VDW(p) = $ ? Is it positive ? 

This question is a "hyperbolic" analog of the van der Waerden conjecture for perma- 
nents/mixed discriminants. If the van der Waerden constant is positive then our analog 
of the Edmonds-Rado theorem for hyperbolic polynomials follows. 
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2.1 "Hyperbolic" scaling 

Let us explain why Question 2 above is indeed an analog of the van der Waerden conjecture 
for hyperbolic polynomials. 

Lemma 2.7: The infimum in (5) is attained iff there exist positive numbers (ai, ...,a n ) with 
rii<j<ri a i = 1 an d an e-positive vector d G C e such that the tuple (a^xi, a n x n ) is d-doubly 
stochastic. 



Proof: The "if part follows directly from Proposition 2.5. Moreover in this case Cap(X.) = 
p(d). The "only if part follows, very similarly to the proof of Proposition 2.5, from the standard 
necessary condition for the corresponding conditional extremum. | 

A fairly direct adaption of Lemmas 3.7, 3.8 from |21| gives that The infimum in (5) is 
attained and unique if a tuple X is e-positive. Therefore Question 2 is equivalent to finding 

inf — p ^ j ' ' — — : (x\, .., x n ) is d — doubly stochastic, d G C e 



Definition 2.8: Consider an e-nonnegative tuple X = (x\, .., x n ) such that the sum of its com- 
ponents S'(X) = d = J2i<i<k x i i s e-positive. Define the following map (Hyperbolic Sinkhorn) 
acting on such tuples: 

^ (x)=Y=( ^k)'-"'^ty ) 

Hyperbolic Sinkhorn Iteration (HSI) is a recursive procedure: 

Xj + i = HS(Xj),j > 0, Xo is an e-nonnegative tuple with X{ G C e . 

l<i<k 

We also define the doubly-stochastic defect of e-nonnegative tuples with e-positive sums as 

DS(X) = Yl ( tr d(xi) ~ I) 2 ; E a* = d € C e 

l<i<k l<i<k 

I 



Example 2.9: Consider the following hyperbolic polynomial in n variables: p(zi, z n ) = 
Yli<i<n z i- ^ is e- hyperbolic for e = (1,1,.. .,1). And N e is a nonnegative orthant, C e is a 
positive orthant. An e-nonnegative tuple X = (xi, .., x n ) can be represented by an n x n matrix 
Ax. with nonnegative entries: the ith column of A is a vector Xi G R n . If Z = (z%, z n ) G i? n 
and d = (di, d n ) G i? n ; z< > 0, 1 < i < n, then tr d (Z) = Ekk„ f • 

Recall that for a square matrix A = {aij : 1 < i,j < A^} row scaling is defined as 

= {^}, 
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column scaling as C(A) = — j. assuming that all denominators are nonzero. The iterative 

process ...CRCR(A) is called Sinkhorn's iterative scaling (SI). In terms of the matrix Ax. the 
map HS(X.) can be realized as follows: 

So, the map HS(X.) is indeed a (rather far-reaching) generalization of Sinkhorn's scaling. Other 
generalizations (not all hyperbolic) can be found in |22| . 0, I 

Before proving the next theorem let us first state and prove some properties of the map HS. 

Lemma 2.10: Consider an e-nonnegative tuple X = (xi,..,x n ) with 5(X) £ C e , i.e. S(X) 
being e-positive. Then Cap(HS(X.)) = (]li<i<n t r S(X)( x i)) C a P@ty> and p(S (H S (X.))) < 
P (S(X)). 

Proof: Consider the following polynomial in one variable D(t) = p(td + x) = X)o<i<n c «^- ^ 
follows from the identity (4) that 

c n = MpCd,..,^!)- 1 =p(d),c n - 1 = M p (x,d,..,d)(V.(n-l)\)-\...,c = M p (x, .., x)(n!)- x =p(x). 

(19) 

Let (X ( f\x) > X^\x) > ... > An^(x)) be the (real) roots of x in the e-positive direction d, i.e. 
the roots of the equation p(td — x) = 0. Notice that 

X ( f\x) = — and tr d (x) = 2llzl M 

l<i<n ° n ° n 

Thus we get the following identities: 

p(x)=p(d) J] X ( f\x),tr d (x) = M p (x,d,..,d)((n-l)lp(d))- 1 (20) 

l<i<n 

The first statement follows directly from the following obvious formula 

Cap{a\xi, ...,a n x n ) = ( J| ai)Cap(xi, ...,x n );ai > . (21) 

l<i<n 

The second identity in (20) reproves the well known fact that the functional tr d {x) is linear. 
Therefore, we get that 

tr s(x)) (HS(X))= £ tr d (-£i-) = n. 

l<i<n % > 

Via the standard arithmetic/geometric means inequality and using the first identity in (20), 
we finally get that 

p(S(HS(X)))=p(d) n \f(S(HS(X)))<p(d)(Y: xf>(S(HS(X))))^=p(d). 

l<i<n l<i<n 

I 

We also need the following "quantitative" version of Proposition 2.5. 
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Proposition 2.11: Consider an e-nonnegative tuple X = (x±, ..,x n ) with d = J2i<i< n x i being 
e-positive. If DS(X) = T,i<i<ki tr d( x i) ~ I) 2 < k then Ca P{^) > 

Proof: Quite naturally in this context (see, for instance, JH]), we will use Proposition 2.2 , 
i.e. we will prove the the conditions of this proposition imply the generalized Edmonds-Rado 
condition : 

Rank(J2i£s x i) — \S\ f° r au & c {!> 2, n}. 

Suppose that the generalized Edmonds-Rado condition does not hold or, wlog, suppose 

(El<i<fcSi) = 0, where 1 < A; < n - 1. 
Since iV e = Nd,C e = C d , we can use Facts 1 and 2 stated in the proof of Proposition 2.2, to 
get that 

Af } ( X! a*) ^ A ? } ( d ) = 1, 1 < » < fc - 1 

l<i<fc 

Therefore, trd(J2i<i<k x i) < k — 1. On the other hand, using the linearity of the functional 
trd(x), we obtain that 

tr d { ^2 %i) = tr d(xi) = ( 1-< ^)> 

l<i<k l<i<k l<i<k 

where Ei<i<n(^) 2 — «■ Therefore, the Cauchy-Schwarz inequality implies that 

£ l*l< 1 /|<i 

l<j<fc v 

This gives that i?*d(X]i<i<fc a?i) > fe — 1, the desired contradiction. 
I 

Theorem 2.12: Consider Hyperbolic Sinkhorn Iteration (HS1) ~X-j+i = HS(X.j); j > 0, 
where the initial e-nonnegative tuple Xo = (xi,..,x n ) satisfies do = J2i<i<n x i 6 C e . T/ien i/ie 
following statements are equivalent: 

1. The exists j > suc/i £/tat DS(Xj) < - 

2. Cap(Xo) > 

5. lim J DS(X i ) = 

Proof: The implication 1 — > 2 is Proposition 2.11, the implication 3 — ^ 1 is obvious. It remains 
to prove that 2^3. Let us introduce the following notations: 

■^-j = \ X j,li ■■■i x j,n), J2l<i<n x j,i = t r dj(Xj t i) = a(j,i), IIl<i<n a {ji V = Fj" 

First, Lemma 2.10 gives that p{dj+\) < p(dj) < p(do). Thus, directly from the definition 
(5), Cap(Xj) < p(do) < oo. In other words, the sequence (Cap(X.j),j > 0) is bounded. 
By the definition of Hyperbolic Sinkhorn Iteration (HSI) we get that 
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Therefore, using (21) and the arithmetic/geometric means inequality, we obtain that 

Cap(X j+1 ) = Ff^apQL&Ff 1 > \-j > 

Moreover, if DS(Xj) does not converge to zero then the product Pj = ]lo<fc<j F~ l converges 
to infinity.Cop(Xj + i) = PjCap(Xo) and Cap(Xo) > 0, therefore if DS(X.j) does not converge 
to zero the sequence (Cap(X.j),j > 0) is not bounded. This is the desired contradiction. | 

Remark 2.13: We can define the map HS(.) directly in terms of the polynomial 
Q(a 1 ,...,a n ) = P Xl ,.., Xn (ai,...,a n ) = p( ^ a»Xi). 

l<i<n 

Indeed, if J2i<i<n a i x i = d £ C e then 

a i ^ r Q(.ai,...,a n ) 

tr d ( ai Xi) = — -i — (22) 

Q{ai,...,a n ) 



As Q{p.\i ...,a n ) is a homogeneous polynomial of degree n thus it satisfies Euler's identity: 

d 

1 da; 



d 

Q(a 1 ,...,a n ) = n 2J oi i —Q{ax 1 ... , a r 

Ki<n 



(Notice that the linearity of tr^ix) for d G C e is a particular case of Euler's identity.) 
Using formula (22), we can redefined the map HS(.) as 

_ Q(ai, .-,««) Q(ai,...,a n ) 
{ax, ...,a n ) — { 9 d i. 

■•■,««) ^(^(ai, ...,a n ) 

Correspondingly, the inequality p(S(HS(K))) < p(5(X)) can be rewritten as 

Q((— Q(ai, ...,a ri ))- 1 , (^—Q(a u a n ))^ x ) < Q(ai, cvT^; a* > (23) 
oai oa n 

where the equality is achieved iff cn-J^- = Q{ol\, ...,a n ). 

The inequality (23) is indeed "hyperbolic" , it is not valid for general homogeneous polynomials 
with nonnegative coefficients. Consider Q(a±, CK2) = af + «§; n = 2. Then 

d d 

Q{(a — Q( a U-> a n))~ l , ■-, (o «n)) ) > Q("l, a n ) _1 

There is another inequality for general homogeneous polynomials with nonnegative coefficients 
involving partial derivatives: the Baum-Snell-Bregman inequality |31|. |32| . |33| . It is interesting 
that in the case of Example 2.9 the Baum-Snell-Bregman inequality is equivalent to (23). Also, 
in the case of Example 2.9 the map HS(.) is a composition of two Bregman's projections 
associated with one convex functional f(xi,...,Xk) = J2i<i<k Xiln(xi); X{ > (pH], |12])- It 
remains to be understood whether the map HS{.) for general hyperbolic polynomials has some 
Bregman's projections interpretation. | 
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3 Obreschkoff theorem and hyperbolic determinant al polyno- 
mials in three variables 



Recall that the companion matrix C q of a monic polynomial q(x) = x n — a\x n 1 — .. — a n is a 
n x n matrix defined as follows: 



( 


1 











1 












\ a n O-n-1 Cln-2 ■ ■ ■ «1 / 

Consider two monic polynomials of degree n, q = x n — a\x n ~ l — ..— a n and r = x n —bix n ~ 1 — ..—b, 
and define the following homogeneous polynomial of degree n in three real variables: 



Definition 3.1: p{x,y,t) = Det(xC q + yC r — tl) | 



Notice that with respect to this polynomial the roots of a vector (x, y, z) with x + y ^ in the 
direction (0,0, 1) are ((x + y)\\ + z, (x + y)X 2 + z, ■ (x + y)\ n + z), where (Ai, A 2 , A n ) are 
the roots of the polynomial xq + yr. 



Proposition 3.2: The polynomial p(x,y,t) is e-hyperbolic, where e = (0,0, 1), iff' all polyno- 
mials of the form {xp + yr : (x, y) / 0; x, y £ R} have only real roots. 

Proof: First, let us prove the "if part. Recall that the eigenvalues of the companion matrix 
C q are exactly the roots of the polynomial q. Therefore, if x + y / then the eigenvalues of 
xC q + yC r are (all real) roots of xp + yr multiplied by a real number x + y. Ifx + y = then 
the eigenvalues of xC q + yC r are (0, 0, xa\ + yb\) ; and thus also real. 

Second, we prove the "only if part. In the view of the first part we need only to prove that 
all roots of polynomial q — r of degree n — 1 are real. Assume, wlog, that q and r don't have 
common roots. Suppose that there exists a complex zq = x + iy, y > such that (q — r){z) = 0. 
In other words the rational nonconstant function R(z) = — 1 has a zero in the upper half- 
plane H + = {z : Im(z) > 0}. Since R{z) is analytic and nonconstant on H + and R(zq) = 0, 
the range {R(z) : \z — zq\ < e} contains a complex ball z : \z\ < 5 > for all small enough e. 
Therefore there exists z\ with Im(zi) > such that = 1 + 6. It follows that the polynomial 
q — (1 + 5)r has a non-real root, but in this case 1 — (1 + 5) = — 5 ^ 0. We got the desired 
contradiction. | 



Corollary 3.3: Consider the following "line" of monic polynomials: P a (x) = aq(x + b + ca) + 
(1 — a)r{x + b + ca), where a G R and b,c are fixed real numbers. Let A a (l) > ... > A a (n) 
be the roots of P a . Let f(x±,...,x n ) be any symmetric and convex on R n functional. Define 
F(a) = /(A a (l), A a (n)). If all polynomials of the form {xp + yr : (x,y) / 0;x,y £ R} have 
only real roots then the function F(.) is convex on R. 
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Proof: In terms of the (0, 0, l)-hyperbolic polynomial p(x, y, t) = Det(xC q + yC T — tl) the roots 
of the polynomial P a are equal to the roots of the vector (a, 1 — a, b + ca) as a + (1 — a) = 1. It 
remains to apply either ^U] or Proposition 1.2. I 



Example 3.4: Consider an arbitrary monic polynomial g of degree n with all real roots, define 
Pa{x) = q(x + b + ca) + aq'(x + b + ca), where g' is a derivative of q. Then r = q + g' is also 
monic, and P a (x) = (l — a)q(x + b + ca)+ar(x + b-\-ca). A well known result gives that that the 
pair (g,r) satisfies the conditions of Corollary 3.3. Let A Q (1) > ... > A a (n) be the roots of P a . 
We get that f k (a) = Ei<i<fc<„ A (i) is a convex function on R, g k (a) = - J2i<i<k<n ^a(n - i) 
is a convex function on R. If q is a monic polynomial q of degree n with n distinct real roots 
and P a (x) = q(x + a) — ag'(x + a), then P a has n disinct real roots for all a £ R. Thus fk(a) 
is differentiable for all k < n ; a direct application of the formula for the derivative of implicit 
functions gives that /{.(0) = 0. Together with the convexity it gives that the global minimum of 
/^(a) is attained at zero, 1 < k < n. It is easy to see that / n (a), which is the sum of the roots 
of P a , is constant on R. Therefore, by a well known result, if ab > and \a\ < \b\ then there 
exists a doubly stochastic matrix D a ^ such that A a = D a ^Ki, (i.e. the vector A a is majorized 



The case of nondistinct roots can be now proved by a standard perturbation argument: if 
a sequence of functions f m : R — ► R, m > 1 converges pointwise on R to a function / : R — > R 
and f m (x) > f m (0);x G R, m > 1 then the inequality fix) > f(0);x £ R also holds. 

The results from this example generalize some results of the recent paper [HOI and simplify 
proofs of others. 

Our solution of Open Problem 3.6 posed in ^U] (see Part 2 of of Corollary 1.3) provides the 
following general majorization result: 

Corollary 3.5: Consider a hyperbolic pair (g, r) of monic polynomials of degree n, i.e. a pair 
(g, r) such that all polynomials of the form {xp + yr : (x,y) ^ 0;x,y £ R} have only real 
roots. Consider two real 3-dimensional vectors X = (x,y,z) and A = (61,62,63). Assume that 
x + y = L 7^ 0, 61 + 62 = M 7^ 0, x + y + <5i + 62 = K 7^ 0. Define the following polynomials 



Let Ax, Ax+Aj be the ordered vectors (from the largest to the smallest) of roots of the degree- 
n polynomials Px , Px+A > -Pa correspondingly. Define an n-dimensional vector ORDx as an 
ordering of the vector LAx, ORDx+a as an ordering of the vector KAx+a, ORD& as an 
ordering of the vector MA^ . Then the vector ORDx+a ~ ORD& is majorized by ORDx ■ 

(Of course, we can now prove many statements of this kind via applying the Second Part of 
Corollary 1.3 in its full generality, i.e. using all Horn's inequalities.) | 



by A 6 ). 



P X (t)=xq(t--)+yr(t 
Px+A(t) = (x + 6 1 )q(t - ^) + (y + 6 2 )r(t 
PA(t) = 6 iq (t-^) + 6 2 r(t- 



z 



L 
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Let us consider two polynomials q, r with real coefficients. Assume that q is monic, the 
degree deg{q) of q is n and also that q has n distinct real roots Ai > ... > X n . If deg{r) < n 
then 

r{z) -A+ £ -^ ]A ,a k eR. (24) 



If a complex number z = u + v i and Re(z) = u, Im{z) = v ^ then 

Call a pair of polynomials (g, r) hyperbolic if all polynomials of the form {xp + yr : (x, y) 7^ 
0;x,y G -R} have only real roots. As explained (and used) in the proof of Proposition 3.2, the 
hyperbolicity of a pair of polynomials (q, r) is equivalent to the property 

im(44) ^ if Jm(z) ^ . (26) 
9(2) 

Therefore if all a k in () are of the same sign, say a k > 0, then the pair of polynomials (q, r) 
is hyperbolic. Moreover in this case n > deg(r) > n — 1 ; if the pair(g, r) is coprime (i.e. no 
common roots) then all polynomials of the form {xp + yr : (x, y) ^ 0;x,y £ R} have distinct 
roots as in this case / = ^jxrj > 1 < A: < n. (Recall that the condition ^7jA > 0, 1 < < n 
forces the interlacing of the roots.) 

What we proved above is a slightly less general (because of the assumption that the roots 
of q are distinct) sufficiency part of the Obreschkoff theorem [33] . We will prove below the 
following analytic version of the necessity part. 

Theorem 3.6: Consider two analytic functions F, G with real Taylor series. Assume that all 
roots of F are real and simple: i.e. the set of roots of F is (X k G R, —00 < L < k < U < 00) 
and F'(X k ) 7^ 0. Assume that all analytic functions of the form {xp + yr : (x, y) 7^ 0; x, y £ R} 
have only real roots. Then > 0. 

G( z) 

Proof: Let H ++ = {z G C : Im(z) > 0} be the upper half-plane. Then -pW is analytic on H ++ . 
Also, as we explained above, the hyperbolicity of the pair (F, G) implies that Im(S4) / 
if z G H ++ . Thus, from connectivity of and continuousness of Im(y^) on H ++ , we 

conclude that ^^(y£jy) has the same sign on H ++ . Say wlog Im(^|y) <0, z£ H ++ . In other 
words — §rr is a Pick function. Therefore it has the following integral representation jSj: 



G(Z) a + bz+ I" l±^-dti(t), z£H ++ (27) 



F(z) 

f(z) " y_oo z 

where a£i?, 0>&£i? and ^ is a nonnegative finite measure on R. Since for all k a real 
number X k is a simple root of F, 

G(X k ) x ,G(A fe + ie) 

hm(A fc + le - X k )- 



F'[X k ) 40 v "^(Afc + ie) 

= innfe f l+%t% t). (28) 
40 J-00 A fc + te-* 
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It is easy to see that 



lim* r 1 + t(Afc + " ) ^) = In-e l+ ^ Xk + ie) d,jt) (29) 



40 J-oo Xk + ie-t e|0 A fc -<J A fc + ie-t 

for all (5 > 0. Using the identity 

1 + tz l + t 2 



z — t z — t 

we get that 



+ t; t,z eC 



limie / — — — -du (t) = hmie / +t)dii(t) 

40 iA fc -5 A fe + ie-t 40 J\ k -s X k + ie-t 



lim!£ / A » +i T i+f_ (M( ) = ib. / Si+i (1+ t ) "!, +< "^- ( Vw 



40 A fc -<S A fe + ie-t 40 J\ k s e 2 + (X k - t) 2 

As fle( ^^^p'^ ) > and the last limit exists and is real, we finally get that > 
0. Notice that if F and G don't have common roots then p7jA > 0. | 

Remark 3.7: It is impossible to have a hyperbolic polynomial p(x, y, z, t) = Det(xC q + yC r + 
zC p — tl) in four variables. Indeed, consider three real monic polynomials q, r,p, all of degree n. 
Then there exists a real vector (x, y, z) ^ such that x + y + z = 01 and deg(Q) < n — 2 , where 
Q = + yr + zp. Assume that q has n distinct real roots (Aj, 1 < i < n). The polynomials q 
and Q have at most n — 2 common (real) roots : 

9 (t) = (x - Ai)...(x - A fc )g(t), Q(t) = (x - Ai)...(x - Afc)Q(t), 

where deg(q) + k = n and deg(Q) + k < n — 2. Therefore , 

l*r = E r n r ;o^a i €/i > * + i<i<n. 

If all / G i?, + 1 < z < n have the same sign, then deg{Q) = n — k — 1 and deg(Q) = n — 1. 
But deg(Q) < n — 2 , therefore / a, £ i?, /c + 1 < i < n don't have the same sign. It 

follows from Theorem 3.5 that there exists z G such that = i G R. This means that 
there exists a linear combination aq + Q with / a 6 I! and some non-real roots. Thus the 
monic polynomial of degree n, P = a _1 ((a + x)q + yr + zp)! has some non-real roots and the 
determinantal polynomial p(x, y, z, t) = Det(xC q + yC r + zC p — tl) is not hyperbolic in the 
direction (0,0,0, 1). I 
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4 More on Alexandrov-Fenchel inequalities for mixed hyper- 
bolic forms 

Let p be an e- hyperbolic polynomial of degree n in m variables. Consider p; x k +i, ...,x n which 
are all e-positive. Associate with them the following polynomial of degree k < n in one variable 

4>k(t) = MJx + tp, ...,x + tp,x k+1 , ...,x n ) 

V v ' 

Then for all x S R the roots of the polynomial </> are all real [23] • Let us call this property 
/c-hyperbolicity. (This essentially follows from the fact that if polynomial in one variable q has 
only real roots then its derivative q' also has only real roots. And the latter fact is a particular 
case of the fact that the (complex) roots of q' belong to the convex hull of the roots of q.) 

The Alexandrov-Fenchel inequalities for mixed hyperbolic forms are directly derived from 
2-hyperbolicity: the discriminant of <f>2 is nonnegative. It is clear (for instance, from ^7j) that 
the Alexandrov-Fenchel inequalities (2-hyperbolicity) are not sufficient to answer Question 1, 
i.e. whether or not 

CO(supp(P xu .., Xn ) n I(n,n) = supp(Q)7 

One possibility would be to use /c-hyperbolicity for all k < n, which gives a lot of other inequal- 
ities |23| expressed in terms of Hurwitz determinants of 4>k and (f)' k . This also might be a way 
to settle Question 2 (van der Waerden conjecture for mixed hyperbolic forms). 

We will finish this section with very direct proof of concavity of ln(p(x)) on the positive 
cone C e . 

Let x,y £ C e and < a < 1. Then 

n\p(ax + (1 — a)y) = —MJax + (1 — a)y, ax + (1 — a)y, ...,ax + (1 — a)y) 
n\ 

= E W^y M ^-r-', (30) 

where M(i) is a mixed hyperbolic form M p (Xj), the n-tuple X contains i copies of x and n — i 
copies of y. It follows from the Alexandrov-Fenchel inequalities that if x, y £ C e then 
< M(i) > yJM(i - l)M(i + 1), 1 < % < n- 1. Thus M(i) > M(0)»M(n) ! ?, which gives the 
following inequality 

M p (ax + (1 - a)y,ax + (1 - a)y,...,ax + (1 - a)y) > (oM(0)» + (1 - a)M(n)n) n . (31) 

Using the concavity of ln(x),x > we get that 

ln(M p (ax + (1 - a)y, ax + (1 - a)y, ax + (1 - a)y)) > aln(M(0)) + (1 - a)ln{M(n)) . (32) 

But M p (ax+(1— a)y, ax+(l— a)y, ax+{l— a)y) = n\p((ax+(l—a)y), M(0) =n\p(x),M(n) = 
n\p(y), so 

ln(p((ax + (1 — o)y)) > aln(p(x)) + (1 — a)ln(p(y)). 



17 



5 Acknowledgments 



First, thanks to the Internet, Google and Pablo Parrilo: I came across hyperbolic polyno- 
mials mainly because of looking for Pablo's e-mail address. And Google provided... the 
rest. For instance, I got the very stimulating preprint [HO] as a result of the Google query 
" +arnold+hyperbolic-|-polynomiar' . 

It is my pleasure to thank Adrian Lewis for numerous as e-mail as well phone communica- 
tions. Thanks to Leiba Rodman for discussions on the subject of Section 3. 

I would like to acknowledge the fantastic library of Los Alamos National Laboratory: all 
references I needed were there. 

References 

[1] J.W. Helton and V. Vinnikov, Linear matrix inequality representation of sets, Tech. report, 
Math. Dept, UCSD, 2002. 

[2] V.Kabanets and R.Impagliazzo, Derandomizing polynomial identity tests means proving 
circuit lower bounds, Electronic Colloq. on Comp. Complex, No. 55, 2002. 

[3] L.Gurvits, Quantum Matching Theory (with new complexity-theoretic, combinatorial and 
topological insights on the nature of the Quantum Entanglement), arXiv.org preprint 
quant-ph/0201022 2, 2002. 

[4] L.Gurvits, Classical deterministic complexity of Edmonds' problem and Quantum En- 
tanglement, arXiv.org preprint quant-ph/0303055 , 2003. (extended version of STOC-2003 
paper). 

[5] L.Gurvits, Classical deterministic complexity of Edmonds' problem and Quantum Entan- 
glement, Proc. 35 ACM Symp. on Theory of Computing, ACM, New York, 2003. 

[6] H.Minc, Permanents, Addison- Wesley, Reading, MA, 1978. 

[7] R. B. Bapat, Mixed discriminants of positive semidefinite matrices, Linear Algebra and its 
Applications 126, 107-124, 1989. 

[8] R. Bhatia, Matrix Analysis, Springer- Verlag, New York, 1997. 

[9] N.V. Krylov, On the genaral notion of fully nonlinear second order elliptic equations, Trans. 
Amer. Math. Soc. (3) 347(1995), 857-895. 

[10] H.H. Bauschke, O. Guler, A.S. Lewis and H.S. Sendov, Hyperbolic polynomials and Convex 
analysis, Canad. J. Math. Vol 53, 2001, 470-488. 

[11] O. Guler, Hyperbolic polynomials and interior point methods for convex programming, 
Math. Oper. Res. (2) 22(1997), 350-377. 



18 



[12] J. Edmonds, System of distinct representatives and linear algebra, Journal of Research of 
the National Bureau of Standards 718, 4(1967), 242-245. 

[13] A. A. Klyachko, Stable bundles, representation theory and Hermitian operators, Selecta 
Math. (N.S.) 3: 419-445, 1998. 

[14] G.P. Egorychev, The solution of van der Waerden's problem for permanents, Advances in 
Math., 42, 299-305, 1981. 

[15] D. I. Falikman, Proof of the van der Waerden's conjecture on the permanent of a doubly 
stochastic matrix, Mat. Zametki 29, 6: 931-938, 957, 1981, (in Russian). 

[16] M. Grotschel, L. Lovasz and A. Schrijver, Geometric Algorithms and Combinatorial 
Optimization, Springer- Verlag, Berlin, 1988. 

[17] L. Gurvits, Van der Waerden Conjecture for Mixed Discriminants, submitted, 2000; ac- 
cepted for publication in Advances in Mathematics, 2001. Also NECI tech. report, 1999. 

[18] L.Garding, An inequality for hyperbolic polynomials, Jour, of Math, and Mech., 8(6): 
957-965, 1959. 

[19] N. Linial, A. Samorodnitsky and A. Wigderson, A deterministic strongly polynomial algo- 
rithm for matrix scaling and approximate permanents, Proc. 30 ACM Symp. on Theory 
of Computing, ACM, New York, 1998. 

[20] L. Gurvits and A. Samorodnitsky, A deterministic polynomial-time algorithm for approxi- 
mating mised discriminant and mixed volume, Proc. 32 ACM Symp. on Theory of Com- 
puting, ACM, New York, 2000. 

[21] L. Gurvits and A. Samorodnitsky, A deterministic algorithm approximating the mixed dis- 
criminant and mixed volume, and a combinatorial corollary, Discrete Comput. Geom. 27: 
531 -550, 2002. 

[22] L. Gurvits and P.Yianilos, The deflation-inflation method for certain semidefinite program- 
ming and maximum determinant completion problems, NECI technical report, 1998. 

[23] A.G. Khovanskii, Analogues of the Aleksandrov-Fenchel inequalities for hyperbolic forms, 
Soviet Math. Dokl. 29(1984), 710-713. 

[24] A.S. Lewis, P.A. Parrilo, M.V. Ramana, The Lax conjecture is true, 2003 preprint. 

[25] E.Pascal, Die Determinanten, Teubner- Verlag, Leipzig, 1900. 

[26] L. Gurvits, Determinantal polynomials, bipartite mixed quantum states, Wick formula and 
generalized permanental inequalities, in preparation, 2002. 

[27] M. Fiedler, Pencils of real symmetric matrices and real algebraic curves, Linear Alg. and 
Appl., 141: 53-60, 1990. 

[28] L. Gurvits and Z.X. Li, Smooth time-periodic feedback solutions for nonholonomic motion 
planning, in Nonholonomic motion planning, ed. by Z.X. Li and J.F. Canny, Kluwer Acad. 
Publ., p. 53-108, 1993. 



19 



[29] L. Gurvits, Control of nonholonomic systems and decomposition of skew symmetric ma- 
trices, in Proc. of CDC-1993, San Antonio,USA, 1993. 

[30] Julius Borcea and Boris Shapiro, Hyperbolic polynomials and spectral order, preprint from 
Shapiro's web page. 

[31] L. E. Baum and G.R. Snell, Growth transformations for functions on manifolds, Pacific J. 
Math. 27, 1968, 211-227. 

[32] L. E. Baum, T. Petrie, G. Soules and N. Weiss, A maximization technique occuring in 
the statistical analysis of probabilistic functions of Markov chains, Ann. Math. Statist. 41, 
1970, 164-171. 

[33] L.M. Bregman, Certain properties of nonnegative matrices and their permanents, Soviet 
Math. Dokl., 14, 1973, 945-949. 

[34] L.M. Bregman, A proof of convergence of the method of G.V. Seleihovkii for a problem 
with transportation-type constraints, Z. Vycisl. Mat. i Mat. Fiz. 7, 1967, 141-156. 

[35] N. Obreschkoff, Verteilung und Berechnung der Nullsttellen reeller Polynome, DVW, 
Berlin, 1963. 



20 



